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Abstract : We propose a model for the geometry of a dynamical spherical shell in which 
the metric is asymptotically Schwarzschild, but deviates from Ricci-flatness in a finite neigh- 
bourhood of the shell. Hence, the geometry corresponds to a 'hairy' black hole, with the hair 
originating on the shell. The metric is regular for an infalling shell, but it bifurcates, lead- 
ing to two disconnected Schwarzschild-like spacetime geometries. The shell is interpreted as 
either collapsing matter or as Hawking radiation, depending on whether or not the shell is 
infalling or outgoing. In this model, the Hawking radiation results from tunnelling between 
the two geometries. Using this model, the back reaction correction from Hawking radiation 
is calculated. 
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1 Introduction 



The back reaction of Hawking radiation on the black hole geometry has attracted con- 
siderable attention during recent years |jl|. However, explicit calculation of the resulting 
corrections to the spectrum have only been performed in the semi-classical approximation 
using a thin shell model for the Hawking radiation field [j^], wherein a "dressed shell" wave 
function was obtained in a reduced phase-space quantization formalism, with the back- 
ground metric fixed to be Ricci-flat, and hence Schwarzschild, everywhere except at the 
location of the shell, assumed to be an imbedded 3-manifold. To lowest order, the shell 
traverses a geodesic in the ambient spacetime. Near the Schwarzschild radius (given in 
terms of the ADM mass Ma by r = 2Ma), the affine parameter t of the geodesic acquires a 
small imaginary part for an outgoing shell. The Hawking radiation is thereby generated in 
a narrow region near the Schwarzschild radius. In fact, the spectrum of Hawking radiation 
is completely determined by Imt, independently of the detailed form of the wave function. 

In this paper we examine the dependence of the correction to the spectrum of Hawking 
radiation on the choice of spacetime metric used to model the back reaction. For this 
purpose, we introduce a more detailed model, so that the minimal modification to the 
Schwarzschild geometry reproduces the expected stress-energy tensor in the asymptotic 
region. The possible choices are severly constrained by this asymptotic condition as well as 
the boundary conditions at the shell. This is reviewed in Section 2. In Section 3, such a 
metric is displayed and its properties are discussed. In the absence of the back reaction, it 
reduces to the Schwarzschild metric, with the shell playing the role of collapsing matter, if 
it is infalling, and of the Hawking radiation, if it is outgoing. In the next approximation, 
the results of Kraus and Wilczek [l| are recovered, but with an additional correction of the 
same order of magnitude. This additional correction arises from the details of the dynamics. 
We also derive the corrections arising from a non-zero shell mass- i.e. from a massive scalar 
field. In Section 4, we examine the stress-energy tensor. In the last section, we conclude our 
discussion and anticipate a more complete quantum mechanical treatment of the problem. 

Our model is based on a reparameterization, r R{r;a{r), (3{r)), where a(f),/3((f) are 
parameters depending on the shell position f. The parameters a(f),/5(f) are uniquely 
determined by the boundary conditions at the shell. The ADM mass is then determined 
from the saddle point condition. The parameter /3(r) satisfies a cubic equation, the physical 
roots of which represent two disconnected spacetime geometries. For an outgoing shell, we 
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find that (3(r) acquires a small imaginary part in the narrow non-classical region around 
the Schwarzschild radius. It is in this region that the Hawking radiation originates from 
tunnelling bwtween the two spacetime geometries. We obtain an expression for Imt from 
the solution of the geodesic equation- the lowest order approximation to the motion of the 
shell. A notable feature of this approach is that the geodesies for an infalling particle, as a 
function of Schwarzschild time t, are complete. This is due to the dynamical distortion of 
the Schwarzschild metric- the back reaction. Hence, the affine parameter of these geodesies 
can be used as the time variable in the calculation of the Hawking spectrum. 

2 Review of Spherical Shells and Hawking Radiation 

We will consider a spherically symmetric geometry in which there is a thin spherical shell 
with mass m. In terms of spherical coordinates (t, r, 0, (p), the metric is of the form 

ds2 ^ _Ar2^t2 ^ ^2 ^ Ar^dt)^ + (dO'^ + sin^ 9 dcj?) , (1) 

where A^, L, N'^ , R are functions of (t, r) only. In fact, the dependence on the time coordinate 
t is only through f(t), where r = f(t) describes the position of the spherical shell. 

The action can be written in ADM form ^ |l| : 



S[f, R, L, N, A^n = J dt pf-MA + J dr (UrR + H^L - NHt - N^'H, 



(2) 



In the above, p^Hr^YIl are the canonical momenta conjugate, respectively, to f,R,L. The 
term M4 is a boundary term which is a constant of the motion and may be identified with 
the ADM mass. We use units in which G = c = 1. The lapse A^ and shift A^*" functions are 
Lagrange multipliers enforcing the constraints: 

Tir ■■= R'UR-LU'L-p6ir-fit))^0. (4) 

These are respectively the Hamiltonian and diffeomorphism constraints. Hatted quantities, 
e.g. L are given by L{t) = L{r,t)\ . Also L = dL/dt,L' = dL/dr, etc. Finally, the 

\r=f(t) 

quantity E := [fp' + m?L'^^ ^ . 

The constraints can be solved for the field momenta [||, ^ 

Hr = Ln'j,/i?', if r/f(t); (5) 
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R [{R'/Lf - 1 + 2Ma/R] ^ if r > r{t) 
R[{R' /Lf -l + 2M/R]^ , \ir<r{t) 



(6) 



Moreover, the discontinuity at the shell gives the following boundary conditions: 



■p/U AR' = -E/R, 



(7) 



where 



AUl := lim (nL(f(t) + e) - UL{f{t) - e)) 



etc. As usual in thin shell formulations, it is assumed that the metric is continuous, but 
the first deriviatives normal to the shell are discontinuous at r = f(t). 

A semi-classical quantization of the system can be performed either in terms of the reduced 



wave function 'ip{t,r) associated with the location of the shell. When expanded in modes 
the wave-function has the form: 



where the index z = 1, 2 is not summed over, and denotes, respectively the quantities defined 
with respect to observers falling with the shell and in the asymptotic region. In particu- 
lar, ai{u)),al{uj) are annihilation and creation operators, while the ui and their complex 
conjugates u* are the mode functions. The mode functions U2, in the asymptotic region, 
are given by U2 = ii(u;, r)e~*'^*. The operators corresponding to the different observers are 
related by Bogoliubov transformations, with the coefficients A{uj,lo'), B(uj,lo') given by the 
Fourier transforms of the wave-function: 



phase space [Q], or via Dirac quantization, wherein the constraints are imposed as operators 
which annihilate physical states |§, ^. In either case, one ends up with a second- quantized 






(10) 



(9) 
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The quantity n{uj) gives the spectrum if the contribution of those particles falhng back 
into the black hole can be ignored. If back reaction is also neglected, n{uj) is given by the 
well-known thermal spectrum js], |6| 

n{uj) = -. (12) 

3 The Back Reaction 



In the absence of the shell, the one-parameter family of Schwarzschild metrics is the unique 
solution of the vacuum Einstein equations with spherical symmetry. The presence of the 
shell disrupts this uniqueness. In fact, classically the metric would be Schwarzschild ev- 
erywhere except at the location of the shell itself. The metric has a jump discontinuity at 
the location of the shell, and the Schwarzschild metrics on either side of the shell can have 
unequal Schwarzschild mass. In a quantum mechanical treatment, classical geometry may 
not be completely established. In a semi-classical analysis, there is a transitory regime of 
duration At specified by the quantum uncertainty principle, during which the Schwarzschild 
mass, or total energy, M is uncertain by AM, where AM At ~ h. We expect that in this 
approximation, the uncertainty in total energy is approximately given by the energy of the 
shell, AM~ na;.|] 

We can model the situation via the "semi-classical Einstein equations" 

Gf,u{x,u;) = 87rT^j.(x,u;), (13) 

where the stress-energy tensor Ti^i^{x,uj) includes the Hawking radiation itself, with a con- 
commitant distortion of the "background" Schwarzschild geometry. The program would 
begin with consideration of some matter field, say a scalar field governed classically by the 
Klein-Gordon equations, propogating in the distorted Schwarzschild background. However, 
in the absence of an exact classical solution which describes the back reaction effect, we 
will propose a physically reasonable metric as a "solution" of Eq.(|l3|). For this purpose, 
we will rely on the constraints and the boundary conditions at the shell to limit the possi- 
bilities. In addition to this, the stress-energy tensor T^y will be required to be consistent 
with well-known predictions for the fiux of Hawking radiation in the asymptotic region, 
i.e., as f — > oo. Moreover, Hawking radiation may be viewed as arising from a mismatch 
between the reference frames of an observer at infinity and that of an infalling observer [^j. 
This mismatch is represented by the Bogolubov coefficients given by Eq. (|To[) . Therefore, in 
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our model, the distorted metric must approach that of an infahing observer as r — > r{t). 
We can satisfy this requirement by choosing coordinates such that in the "inner region" 
r < r{t), n^, = = II/j. These conditions severly Umit the choice of metric. 

We have found the fohowing metric satisfies the above requirements: 

f :=r[l+j(r,f(t),uj)], if r > f(t) 
R=\ (14) 
[ R_ ■.= r[l + (3{f{t),uj)], if r<f(t) 

^^i L+:={1-2Ma/R+)-K ^ if r>f(t) 

\ L_ :=(l + /3(f(t),w))(l-2M/i2_)-5, ifr<f(t) 

Furthermore, we wih choose coordinates such that the lapse and shift functions are = 

1/L, iV^ = 0. 

The functional dependence of the parameters j{r,f,uj),(3{f,Lo) on f will be determined by 
the boundary conditions at the shell and in the asymptotic region. 

In the asymptotic region, given by r ^ oo, we require that R^ r. Hence 

7(r,f,u;)^0, (16) 

in the above limit. This ensures that the metric is asymptotically Schwarzschild. 

The continuity of R at the shell requires that ^{f^f^uj) = P{f,Lo). Furthermore, the deriva- 
tive of j{r,f,uj) with respect to its first argument, evaluated at r = r is determined from 
the junction conditions Eq.(|^ and the constraints Eq.(^),Eq.(^). The result is: 

™ = i^ + (17) 

2(l+/3) ^ 2/3(1 +/3-2M/r)' ^ ^ 

where we have written: 

-1 d 

p(r, Lo) or \r=r 

We specify the function 7(r, f,a;) only up to the correct limiting behaviour, as discussed 
above. 

It remains to ensure that L is continuous across r = f{t). This will determine the parameter 
(3. The equation = L_ at r = r{t) implies that 

- {hf)x^ - X + / = 0, (19) 

where 

x:=l + P; f := 2M/f; h := Ma/M. (20) 
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The three roots of this cubic are given by 



hf P+ + P- ^fP+-P\ hf , , 



where 



a= -l-\{hff- h = -^{hff-\[hf)+f- (22) 
V4+27= 

P.= (-^±c) . (24) 

If > 0, there is one real and a complex conjugate pair of roots; if < 0, there are three 
real roots. 

The nature of the roots depends on the value of h. In fact, /i = 1 it tu/M, with the (+) 
sign for ingoing, and the (— ) sign for outgoing radiation. It is easy to see that all the 
roots are real for ingoing radiation (h > 1), but complex roots are possible for outgoing 
radiation. Figures la and lb display the behaviour of the roots, for infalling and outgoing 
shells, respectively. ||7| 

The behaviour of the roots in the respective regions where f » 2M or where f ~ 2M is 
most important. In the former, where the shell is far from the horizon, x is very close to 1. 
We call this root "physical", since for it, the metric is asymptotically flat. We note that for 
the physical root, /3 = ±uj/f + O (r~^) for the + (— ) sign for ingoing (outgoing) shell. 

In the other region, near the event horizon, we expand in powers of the two parameters 
uj/M and e := r — 2M. To lowest order, for the outgoing shell, the roots are 

where 



A:: ^ 
2M 



" M; 2 \2M 



(26) 



To summarize, the physical root 1+/3 of Eq.([l9|) is always real for an ingoing shell, but for an 
outgoing shell, it acquires a small imaginary part in a narrow region around r = 2M. This is 
the region where Hawking radiation originates and escapes by way of tunnelling through the 
classical barrier. In this region the solution of the geodesic equation can become complex, 
as will be shown below. 
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Figure la: x vs. Shell Position f,h= 1.01 



At the classical level, each point on the spherical shell traverses a geodesic. The geodesic 
equation can be obtained directly from the boundary condition Eq.(|^ and from the expres- 
sion for the momentum canonically conjugate to r, 



p = mL r 



(27) 



The result is that 



f = ±N (aUl/R) /AR' 



(28) 



where the + (— ) sign applies to an outgoing (ingoing) geodesic. (If m = 0, it follows 
from Eq.(28) and from Eq.(|^) that f = =biV, so that the points on the shell traverse null 
geodesies.) 



The geodesic equation Eq.(|28|) is integrated to yield 



dfR 



^-2^^2(1 + W(i?;)^-1 



R ( m 



R-2M \r 



+ to, (29) 



where R = f{l + P) and rg =constant. For /3 = and m = 0, this reduces to the usual 
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Figure lb: x vs. Shell Position f,h = 0.99 



equation for a null geodesic in Schwarzschild spacetime: 

t = ±[f + 2M\n{f/2M-l)]+to. (30) 

However, if f3 is complex, then t becomes complex. In this case, the integral has a branch 
cut along the real axis in the complex f-plane. The region in which f3 is complex is quite 
narrow, centered around f = 2Ma- In this region, /3 « constant, and the branch cut can be 
replaced with a simple pole. The imaginary part of t can then be evaluated as the residue 
at the pole, with the result: 



^ 16 \MaJ 



(31) 

Here we have also assumed that m « lo. Only the lowest order corrections in lo / Ma, m/M a 



1 / m \2 



and m/(jj are retained Eq.(31). If m = 0, one obtains 



Imt = 47rMA(l - 3uj/8Ma) = A'k{Ma - 3u;/8). (32) 

If the last term is ignored, we have the result of Ref. 1: Imt = AttMa = 4:Tt{M — lo). If all 
the corrections of order uj/M are ignored, we have Hawking's result: Imt = 47rM, obtained 
in the Euclidean formulation. ||5|, |^ 



The spectrum of Hawking radiation is given by Eq.(|T^) and Eg. ([Tl|) . By virtue of Eq.(|3l[) 
we obtain: 

f r I Tn2 \ 1 1 

■ - 1. (33) 



[n(a;)] ^ = exp |87ru; 
For the case of a massless shell (m = 0), we have: 



3 / , m 

Ma - 1 + 



16M| 



= exp [Svrcj (Ma - 3w/8)] - 1. 



(34) 



With Ma = M — Lij, we recover the results of Ref. 1, if the additional correction of 3cj/8 
is ignored. The origin of this additional correction lies in the details of the dynamics, 



which are sensitive to the form of the metric. Our choice of metric, given by Eq.(14) and 
Eq.(^), reproduces the results of Ref. 1, if the details of the shell dynamics in the distorted 
background, represented by Eq.(^) are ignored. However, this does not seem to be justified 
since the additional correction in Eq.(p3|) is of the same order in a; as the corrections in 
Ref. 1. It appears that the back reaction correction is sensitive to the details of the metric. 
The ambiguity in the metric can only be diminished by an exact semi-classical solution of 
the Einstein equations. In the absence of such solutions, we have to rely on the consistency 
of our choice of metric. We will confirm this consistency through a calculation of the 
stress-energy tensor Tp^y in the next section. 



4 The Stress-Energy Tensor 

We can calculate T^^y using the standard thin shell formalism Q . We use Gaussian normal 
coordinates (r/, x*, i = 1, 2, 3) so that the metric is of the form: 

ds^ = drf' + hijdx^dx^ , (35) 

where r]{x'^) = is the equation of the shell and hij is the induced metric on the hypersurface 
Sj. swept out by the trajectory of the shell. In our case, rj = r — r{t) for a spherical shell. 
The stress-energy tensor takes the form 

= siv) + 0(7?) r+ + ei-7j) (36) 

where 9 is the unit step- function. The tensor S^j^u is the surface stress-energy carried by 
the shell, while denote the regular background contributions on the two sides of the 
shell. Hence, describe the black hole hair originating on the shell. Both fall off rapidly 
with distance from the shell. Therefore Hawking radiation contributes to Sfj,^, but not 
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directly to T^. The latter represent the transient stress-energy originating in the quantum 
fluctuations of the geometry. 

The surface stress-energy tensor S^y is determined by the extrinsic curvature (second fun- 
damental form) K^jy of by 

where is the unit normal to Tir- The Einstein equations now take the form: 

Sij = - (l/Svr) {Kij - Khij) , (38) 
DjS'^ = -r;, (39) 
^[k+ + k^)s^^ = Tl (40) 

where Kij := Kf^ — K^j, k := h''^ Kij and Dj is the covariant derivative with respect to the 
metric hij. 

We are primarily interested in the component Su for an outgoing shell, which determines 
the surface energy per unit area of the shell. Other components of the surface stress-energy 
tensor determine the pressure and tension. Prom Eq.(37)-Eq.(|38[) we obtain for a massless 
shell: 

Su = -N^Ar'/AttR, (41) 

(42) 



as f — > oo. Therefore, energy outflow of Hawking radiation from a sphere of radius r at 
infinity is huj, as expected. As a result, Eq.(|l4[)-Eq.(^) are consistent with the expected 
asymptotic behaviour of the stress-energy tensor. Moreover, it can be shown that the 
fall off rapidly away from and therefore do not contribute to the Hawking radiation. 



The total flux of Hawking radiation can now be calculated using Eq.(12). The total flux is 
/(^(27r)~^(iw n(u;)Lj = l/7687rM^, which recovers a well-known result Q. 

5 Analysis of the Geometry 

As we have seen, the spacetime geometry is determined by the properties of the function 
P{f,ui). This function must satisfy the cubic equation Eq.([l9|). The properties of these 
roots are displayed in the graphs of the 'optical scalar' R vs. f. Figure 2a and Figure 2b. 
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Figure 2a: Optical Scalar R vs. Shell Position f,h = 1.01 



We may conclude from these that the root for which < for all f > is unphysical, while 
the other two roots of the cubic represent a bifurcation of the spacetime geometry. 

The collapsing shell, which has Ma = M + for which the optical scalar corresponds to 
a Schwarzschild geometry asymptotically, i.e. which behaves as i? r for large r, never 
reaches the location R = 2Ma except in the limit as r ^ 0. The latter is a singular event 
in the geometry. The events at ii = 2Ma and at -R = 2M constitute apparent horizons. 
The two geometries are separated by a gap 2M < R < 2Ma- A collapsing shell inside the 
R = 2M horizon does fall into the singularity in a finite time. However, it cannot cross the 
apparent horizon R = 2M in a finite time. Classically these two configurations are distinct 
solutions, and represent different spacetime geometries. It appears that the shell will only 
collapse into a black hole via quantum fluctuations of the spacetime geometry. Similarly, a 
shell in the inner region geometry can enter the outer region only via quantum tunnelling. 
Since the magnitude of the distortion of the outer geometry from that of Schwarzschild is 
small in the region where f is near the value 2Ma', and since, moreover, t is continuous at 
r = 2Ma, t can be used as an afiine parameter. 
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Figure 2b: Optical Scalar R vs. Shell Position f,h = 0.99 

When Myi = M — lo, i.e., the shell is outgoing and the metric is complex in a small region 
near where the shell is close to the horizon.]^] A classical spacetime geometry does not exist 
in this region, and the Hawking radiation from classically forbidden configurations takes 
place through quantum tunnelling. 

To summarize, in our model the fixed Schwarzschild geometry is split by the back reaction 
into two distinct spacetime geometries, both of which contain apparent horizons. The 
exterior geometry approaches the schwarzschild limit as r ^ oo, while the interior geometry 
is Schwarzschildean in the limit f ^ 0. 

6 Conclusion 

We have calculated the back reaction correction for a massless as well as a massive thin shell 
in the semi-classical approximation. We have found an additional correction due to the back 
reaction for a massless shell, so that the total correction to M in the Hawking spectrum is 
— lla;/8 as opposed to —uj, obtained in Ref. 1. The correction appears to be sensitive to 
the details of the dynamics, and, therefore, to the choice of metric distorted from that of 
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Schwarzschild by the presence of the sheh. Our choice satisfies the asymptotic conditions 
as weU as other consistencey conditions at the shell. However, the ambiguity in the choice 
of metric will remain until an exact solution of the semi-classical Einstein equations can be 
found. 1^] 
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